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Abstract 

A microscopic description of the quasielastic {p,n) charge-exchange reaction (here, charge- 
exchange scattering between analogue states) is presented and discussed. Emphasis is focused 
on the spin-isospin structure of the projectile-target coupling. The model is a coupled-channel ex- 
tension of the full-folding optical model approach (OMP) developed for nucleon elastic scattering, 
where emphasis is placed on retaining the genuine off-shell behavior of realistic effective interactions 
in the nuclear medium. The resulting non-local optical potentials are applied to the calculation 



^3 I of {p,n) differential cross sections, with particular emphasis on small-angle Fermi (A^* = 0) cross- 

■ sections to isobaric analog states. These parameter-free results provide a reasonable description of 

the ^^C(p,n)-data at proton energies above ~100 MeV, but deteriorate for heavier targets. These 
^ ' shortcomings are analyzed and possible ways to correct them are discussed. 
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I. INTRODUCTION 



An issue of high current interest in nuclear research is the physics behind systems and 
processes where the isovector component of the nuclear interaction is relevant. This par- 
ticular aspect is of pivotal importance in understanding a wide variety of phenomena such 
as extreme isospin asymmetric nuclear matter, nuclear systems far from stability, collisions 
of radioactive beams, the asymmetry energy of nuclear systems and {p,n) charge-exchange 
processes A sound understanding of these phenomena may have signifi- 

cant impact on other areas of research such the physics of neutron stars, the formation of 
astrophysical objects, etc. Our focus here is on the study of charge-exchange reactions at in- 
termediate energies. This effort attempts to provide a rigorous microscopic approach for the 
study of these reactions with emphasis on a detailed account of genuine energy, momentum 
and density effects in the nucleon-nucleon (A'^A'^) effective interaction. 

The study of {p,n) charge-exchange reactions at intermediate energies has received con- 
siderable attention during the past few decades. These reactions have been of great value in 
understanding the isovector modes of excitation of nuclei. At beam energies above 100 MeV, 
nucleon charge-exchange reactions can be considered as a one step process, thus allowing a 
rather clean separation of the nuclear structure from the underlying NN effective interac- 
tion. These arguments have resulted in the usual tp structure of the coupling to the nucleus. 
Furthermore, the weak strength of the isovector component of the interaction relative to 
its isoscalar counterpart has made suitable the use of the distorted wave Born approxima- 
tion (DWBA). In spite of these advantageous considerations, no microscopic effort has been 
able to satisfactorily describe the differential cross-section data without phenomenological 
adjustments a posteriori of the nucleon-nucleus {NA) coupling. This impediment is quite 
significant in that it points to a lack of understanding of the simplest process beyond elastic 
collisions. Therefore, a comprehensive and consistent re- analysis of the problem is relevant. 

Studies of nucleon elastic scattering at intermediate energies (i.e. between 100 MeV 
and 1 GeV) have demonstrated the importance of off-shell effects as accounted for within 
the full-folding OMP for (NA) scattering. In particular, it was shown that the tp scheme 
(as most often applied) is insufficient to satisfactorily describe intermediate energy data. 
Furthermore, nuclear medium effects as included by using nuclear matter g matrices have 
provided significant improvements in the description of the elastic scattering data aaa. 
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The issue which naturally arises is whether these improvements within the full-folding OMP 
to NA elastic scattering are also important in charge-exchange processes. 

Although significant advances have been made with the introduction of full-folding op- 
tical potentials in the description of nucleon elastic scattering, none of these improvements 
have been included in charge-exchange reactions |[10]. Furthermore, most of the reported 
{p,n) reaction calculations, which provide a reasonable description of the data, are semi- 
phenomenological in the sense that well fitted potentials in the elastic channel are used 
to obtain the charge-exchange cross-section by means of the DWBA jlll . From this 
prospective additional theoretical effort is needed to clarify the issues which prevent us from 
obtaining a unified understanding of these phenomena. Along this line, the aim of the 
present work is to asses the importance of treating the intrinsic off-shell behavior of, and 
medium corrections to, realistic nuclear effective interactions in the calculation of nuclear 
charge-exchange reactions. For this purpose we present the spin-isospin formalism needed to 
calculate these couplings within the non-local full-folding optical model and have developed 
the required coupled-channel scattering codes to solve the scattering problem exactly. Our 
work is focused on the AS = isobaric analog state (IAS) transitions for various targets 
and incident proton energies. 

We organize this article in five sections as follows. In Section II we review the spin-isospin 
structure of the NA coupling to display its formal structure suited for charge-exchange 
reactions from a finite nucleus. In Section III we outline the coupled-channel approach for 
studying the {p,n) charge-exchange reactions with non-local potentials in the presence of 
the Coulomb interaction. In Section IV we present and interpret results for representative 
applications and in Section V we summarize the present work and draw the main conclusions. 

II. OPTICAL POTENTIAL: SPIN AND ISOSPIN CONSIDERATIONS 

For simplicity in this discussion, we use the free t matrix as the NN effective interaction 
and consider nonrelativistic kinematics. Both of these restrictions will be removed in the 
following sections. Thus, the generalized optical potential for nucleon scattering and charge- 



3 



exchange reactions may be written in momentum space as 

(kV V I U I kuf,) ^ II '^P'^^P^^ I ^ln.',in'(P') V^in.,in(P) I I) 



m.nv ,n,n' 

I / „/ \J , ,1 ..I 



{p'm' n\\s! p' iJ \t I pmn, ki//^)^. (1) 

Here {t) ^ is the antisymmetrized NN t matrix; u and /i denote the initial spin and isospin 
projections of the projectile and 

where ■0i-m,i-n(p) annihilates a nucleon with momentum p and spin and isospin projections 
— m and —n, respectively. The choice of the pair {/i, /i') is determined by the reaction being 
considered. In the convention in which the proton has isospin projection +i, 

(a*) A* ) ~ d' 2)' d' "i)' (~2' 2)' (~2' "i) 

for the {p,p), {p,n), {n,p), and {n,n) reactions, respectively. Following some recoupling to 
display the transferred quanta and the NN spin and isospin, U becomes: 

Sp, ,Tp 

II dp dp'{F I [V't i(p') ® V^ii(p)]°°^'^^"* I I) 
J2 {SSpM SM'){-y-^{2T+l) 

S,T, M, M' 

W{im;TpT){p'k',SM',T\t \ pk,SM,T)A, (2) 

where 

[V'ii(p') ® V^U(P)]''''^™ = E (V> I SM)nn'ln I TAT) V'L, . ,(p') ^MniP) , 
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2 2 



(3) 

defines the order of coupling. Tp and Tf denote the isospin transfers to the projectile and 
target, respectively. Similarly, Sp and St denote the corresponding spin transfers. Also, 

§(0,0) = 1, S(l,-i/p) =a_,^, T(0,0) = 1, T{1, -nt)=T.n,, 

where a^^,^ (T-nt) is the —Up {—rit) spherical component of the Pauli spin (isospin) matrix. 
The sum over M, M' projects out that part of the t matrix which acts in the spin state S 



and which is of rank Sp in spin space. In particular, we write the t matrix as a sum over the 
contributing ranks (A) in spin space 

t = tj:i'T-Of^Vr, (4) 

A=0 S,T 

where 

0^^ = Vs ho + (<?! + ^2) hi hi + (<?i ® 0^2)^ h2 ■ (5) 

Here, Vs and Vt are the projection operators onto the NN states of spin S and isospin T, 
respectively. The optical potential U corresponding to zero spin transfer to the target and 
Tt = Tp isospin transfer to the target and projectile becomes: 

{u' I S{Sp,u' -u) I u) 



2Sp + l 



(-)i-^(2T+l)(2^ + l)H^(iiii;T,T) (p'k' | 4'/ \ pk)^ . (6) 

ST 

The nuclear structure is contained in the mixed transition density defined by 

P^V(P',P) = 2 {F I [^li(p') ® V^ii(p)]°°^'^' I /) , (7) 

22 22 

which is an isoscalar or isovector as Tp = is or 1. With this normalization for p^j , 
J dp p5/(P)P) = ^ ) the number of nucleons in state I. In the evaluation of the optical 
potential it is useful to make explicit use of the translational invariance characteristic of the 
free NN interaction (or that in infinite nuclear matter) by writing it as 

(p' k' I 4^.'') I p k)^ = 5(p + k - p' - k') {k,' I tg'^ I k)a , (8) 

where the relative momenta k and /t'are given by: 

P — k , p' — k' 
K = , k' = . (9) 

For brevity, it is also convenient to introduce the participating part of the NN interaction 
from Eq. ([6]) as 

{k' I t^^^' I k)a = i^(_)i+^(2T + 1)(25+ 1) T,T)(a.' I tg^^ I , (10) 

S,T 

where W is the usual Racah coefficient, 
and St = for the optical potential. 




A. Scattering from a Fermi Gas (iV > Z): an example 



For elastic nucleon scattering, the transition density of Eq. (jTj) becomes 

p^V(P , P) = 2 5(p' - p) [eiPFin) -p) + i- f^ Q{pp{u) - p)] , (11) 

where G(x) is the Heaviside step function and Pf{t^) and Pf{v) denote the Fermi momenta 
for protons and neutrons, respectively. The transition density appropriate for the excitation 
of the analogue state via the charge-exchange reaction is 



PfV(p ,p) = ^Ij^ '^(P - P) [Q{PF{y) -p) - Q{pf{^)-p)] , (12) 



where | F) satisfies 



' -Z ' ' 



In the Fermi gas (FG) approximation we must have from Eqs. ( l8|9|TTi) . p = p', k = k' and 
therefore n = k' . This insures that only the Sp = part of the t matrix participates as the 
spin-orbit and tensor terms do not contribute to forward scattering in the NN system. With 
these observations, the optical potential for elastic scattering given by Eq. ([6]) becomes 

Ufg = 2 5,,,, (5(k- k') ^ rfp {0 [pf{ /) - p ] (ac | t° + ^^)a + O [pf{u) - p) ] (/^ | t° - i^)a] 

(13) 

where / {u) denotes target nucleons like (unlike) the projectile; e.g. for proton scattering, / 
corresponds to target protons and u to target neutrons. Also, we have dropped the Sp = 
superscript on the reduced t matrix t"^' . Separating the isoscalar and isovector contributions 
to Eq. ([T3]) gives 

Ufg = 2 S^^^> 5^^^> 5(k- k')x 

j d^{{Q[pF{l)-p\ + Q[pF{u)-p\}{t^ I t°|Av)A 

+{Q[Pf{1)-p] - Q[pf{u)-p\}{k I I . (14) 

For the nucleon charge-exchange reaction exciting the isobaric analogue state, the Fermi gas 
transition potential is 



-4 ^ 

Ufg = ^ ^1^,^' ^(k - k') 

j dp {Q[pF{y) -p] - e[pp(7r) I I h)a. (15) 



Some insight into the relationship of the optical model in the Fermi gas to that in finite 
nuclei may be obtained by using the Fourier representation of the momentum-conserving 
delta function appearing as a factor in Upc above, e.g. 



5(k-k') = ^/rfre'(-'^')-, 



together with the assumption that the t matrix varies little over the range of p allowed by 
the Fermi momentum. When this is done, the integrals over the step functions may be done 
giving 

for protons (vr) and neutrons (z/), and the optical potential for elastic scattering becomes 

Ufg = S,,^' j dv e^(k-k') r ( I tO + ti I n)A + Pu (a^ I - tM k)a) , (16) 

where n is to be evaluated at p = , for example. In the strict Fermi gas model only the 
exponential above depends on r and we recover the delta function 5(k — k'). When k = k' 
the integration over r gives the volume VL occupied by the gas and the product VLpi gives the 
number of target nucleons of type i. Therefore, 

Ufg{\^ = k') = 5^,^, 5,,,, ( iV;(A^ I t° + i^)a + iV,(A^ I t° - a^)a ) , (17) 

where Ni {Nu) is the number of target nucleons like (unlike) the projectile. Similarly, the 
non- vanishing diagonal elements of the optical potential for charge-exchange reactions is 



Ufg{\^ = k') = -2VN^ ^u,u' (a^ I i^)a . (18) 

If in Eq. ( IT6l) . we use the local density approximation for by allowing it to be a function 
of r, then 

Ufg = 5;.,^' 5.,.' ( Pi{q) (k I t° + t^)A + Pu{q) | t° - | k,)a ) , (19) 

where q = k — k', the momentum transferred to the target and pi is the Fourier transform of 
the density of nucleons of type i. This is an approximate form of the usual tp approximation 
to the optical potential. 
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B. Scattering From Finite Nuclei 



Results for finite nuclei may be obtained by expanding the creation and annihilation 
operators of Eq. ([3]) in a (spherical) shell-model basis. The result for arbitrary (allowed) 
spin and isospin transfer is 



^ii(p')®V'u(P) 



SM;TN 



^ Unljip) Un'l'j'ip') [L Jjj' 

nlj,n'l'j' L,J 



\l'+L+S+J] 



(v \ f\ 



L S J 



{[r,(p')®rKp) 



J-TN iSM 



(20) 



where, in this equation, L, S, J, T denote transferred quanta Lt, Tt etc. Also, c'^ and c denote 
creation operators for particles and holes, respectively, and Unij{p) is the radial part of a 
shell-model orbital in momentum space. The symbol in rounded parentheses is a nine-j 



symbol and ct denotes \/2x + 1. Inserting Eq. fl20l) into Eq. ([7]) gives the transition density 
for a finite nucleus 



2 ^ Unljip) Un'l'j'ip') X] [^tJtJj'i-Y' 
nlj,n'l'j' Lt Jt 



+Lt+St+Jt 



J 

Lt St Jt 



StMt 



(21) 



Although only the 5'f = part of the above density is required for the usual case of calculating 
the optical potential, it is useful to have the more general result for later use. For the special 
case in which the transferred quanta iLt St Jt) are either assumed or restricted to be zero, 
Eq. fl2Tl) becomes 

-V2 



P?/(P',P) 



47r 



E 



tlVV 



^ Unljip) Unljip) Plip' ■ P) 



nlj 



n'lj-n;^u' nlj-fi; 



I) . (22) 



For the special case of elastic scattering in which \ F)=\I),n = n' and N = , the 
mixed density reduces to 



PfI ^(p'' P) = ^ X] ^riljip) Unljip) [Znlj + (-)^* Nnlj\ P/(p " p) 
nlj 



(23) 
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where Znij and Nnij denote the number of protons and neutrons in the state (n/j) where it 
has been assumed that the orbitals for neutrons and protons are the same. For the {ji,n) 
charge-exchange reaction to the isobaric analogue state, z/ = i/' = |, Tt = = 1 and the 
Clebsch-Gordan coefficient is unity giving 



Plr\p'^ P) = - V E ""-M UuiAp) [Nnl, - Znl,] m' ■ P)- (24) 

nlj 

From Eq. ([6]), the optical potential for elastic scattering associated with the transition 
density of Eq. (!23|l reduces to 



' pp,nn 



I LI ij unijyf) unijyf ) 1 lyiJ ■ iJ) ^ ^ 

nlj 



2Sp + l 



(AT^ (k' I t^^^ + t^^^ I k,)a + iV^. {k,' I t^-° - t^^'i I Ai)^) , (25) 

where iVj;^- and A^^^- denote the number of target nucleons in the state (nlj) which are like 
and unlike (isospin projection) the projectile and p' = p + q. Similarly, the optical potential 
mediating the {p,n) reaction to the isobaric analogue of the target ground state is 

2 



tl-l,^l' 



nil 



^ j dp Unljip) Unljip') i ^nlj " Znlj) Plip' ■ p) 



The above expressions can be further simplified in the case of non-spin transfer to the 
projectile. Denoting the proton and neutron mixed densities by Pp and p„, 

Pp(p',P) = Unljip) Unljip) Znlj Pliff' ■ p) , 

nlj 

Pn(p',P) = 5Z ""'i(P) ^nljip') Nnlj PliP' ' p) , (27) 
nlj 

we are left with the following three terms to evaluate 

Upp = ^ j dp [Ppip', p) {i^' I t° + I k) A + Pn(p', p) I t° - I k)a\ , 

Ux = - — 7== dp [pnip', p) - ppip', p)] ii^' I I «;) A , 
/vrv ly — Z J 

Unn = ^jdp[ppip',p){K,'\t'-t'\K,)A + Pnip',p){l^\t' + t'\H,)A\ • (28) 
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C. In-medium calculations 



The optical potential Upp, Unn and Ux obtained in the previous section requires the 
convolution of single-particle wavefunctions with a two-body effective interaction. Simpler 
expressions for these potentials are obtained with the use of the Slater approximation to the 
ground-state mixed density in the case of spin-saturated targets. This approximation was 
examined in Ref. [l3| in the context of nucleon elastic scattering at intermediate energies, 
and its effects are reported to be only noticeable at the larger scattering angles. Furthermore, 
this representation of the mixed density yields a simple prescription for including explicit 
medium effects by means of infinite nuclear matter g matrices. Following the same argu- 
ments discussed in Ref. {q] to incorporate nuclear medium effects, the first of the potentials 
expressed in Eqs. f l28|) takes the general form (c.f. Eq. (21) of Ref . [3]) 

Upp{k', k) = jdR e^('»'-'=)-« [pp(i^) g^^{k', k) + p^{R) g^p{k', k)] , (29) 

where pp and pn are the local proton and neutron point densities, respectively, and (Jnn rep- 
resent off-shell Fermi- averaged amplitudes in the appropriate NN channel. More explicitly, 
in the context of infinite nuclear matter these density- dependent amplitudes are given by 

g^^ik', k) = ^ [eCk- |P|)(?k+p(a^', K;V^;p)dP, (30) 
47rA;3 J 

where gK+p{i^' ■, i^', y/s] p) is taken as the g matrix for symmetric nuclear matter of density 
p. Here, 

K=\{k + k'), P = l{p + p'), k + p=k' + p', (31) 
and the relativistically corrected relative momenta k and k,' are given by 

K = Wk-{l- W)p , k' = W'k' -{1- W) p , (32) 
with W and W scalar functions of the momenta of the colliding particles with relativistic 

n . . 

kinematics built in |7|. Additionally, the local momentum k is taken from A;^ = 37r^p/2, with 
p [Pp{R) + Pn(-R)]/2, the isoscalar local density at radius R. Analogous expressions are 
obtained for Ux and Unn- 



III. COUPLED-CHANNEL CALCULATIONS 



Our study focuses on quasielastic scattering to the isobaric analog state. Considering 
explicitly the isospin degrees of freedom of the scattering waves in the form of outgoing 
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jroton and neutron wavefunctions, we obtain a non-local version of the coupled- channel 



14 



151 ] equations 



Kp + Uj,i + Vc 



(33) 



Here and f/^ are the non-local potentials given by Eqs. (!25|) and (!26|) in the previous Sec- 
tion. Furthermore, Upp represents the short-range coupling between the charged projectile 
and the target protons, i.e. the hadronic-plus-Coulomb contribution with the point- Coulomb 
potential, Vc = Ze^ /r, subtracted. A formal solution to these equations is given in the form 
of the Lippmann-Schwinger integral equation with proton waves $p in the entrance channel 







GpUpp GpUx 

where Qp and Qn are the Green's functions for outgoing protons and neutrons, respectively. 
In coordinate space, for partial wave /, these propagators are expressed in terms of the 
Coulomb (and Bessel) spherical waves 

i 2e 



(34) 



Gl^\r,r';k) = -y ^ Fi{r]; kr^) [Fi{ri; kr^) + iGi{ri; kr^)] 
k 

= -j:^ Fi{r];kr^)Hj^\r];kr^^ 



(35) 



Here 1] = Ze'^/hv, represents the charge parameter for the NA coupling, with v the corre- 
sponding relative velocity, and e the ejectile-residual nucleus reduced energy. For uncharged 
particles Fi and Gi become the usual Ricatti-Bessel functions with the following the phase 
conventions: Fi{0;t) = tji(t); jo(t) = sint/t, and G/(0;t) = tni(t); no(t) = — cost/t. 

The primary input in these equations is the optical potential which we obtain following 
a current version of the full-folding OMP to nucleon scattering . In this scheme, the 
optical potential is calculated in momentum space leading to a general non-local structure 
in coordinate space, a feature which is retained throughout. An exact treatment of the 
Coulomb interaction in the presence of this non-local coupling is performed in coordinate 
space by solving Eq. flMl) once the Fourier transform of the hadronic contribution have been 
performed. Thus, the scattering amplitude is readily obtained from the asymptotic behavior 
of the radial scattering waves from Eq. which, schematically, has the form 





















Un,l 


(r)_ 





^Fi{ri;kpr) 



+ 



^PpiHl^\v;kpr) 
A,„^ifW(0;A:„r) 



(36) 
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with kp and kn the momenta of the outgoing protons and neutrons, respectively. The above 
result provides alternative forms for obtaining the scattering amplitude, namely by simple 
identification of the asymptotic form of the solution, or by explicit evaluation of the matrix 
elements, i.e. 



^pp ( ^pp ^Pi^ ~^ ^n,l ) — ikp fpp , 

A('^ = ( Up,i + Unn Un,l ) = ikn f^l . (37) 



More explicitly. 



{Uifupj) = -t^ dr dr'Fi{7];kpr')[r'Uj,}{r',r)r]up,i{r). (38) 

Analogous expressions hold for ( Upn Un,i ) and ( Unn Un,i ) ■ 

The results obtained with this procedure were compared with DWBA results and showed 
almost unnoticeable differences. 



IV. APPLICATIONS 

The scattering calculations performed in this work differ from most reported applications 
of quasielastic nucleon exchange processes to date. Thus, it is important to spell out the 
differences explicitly. The optical potential is obtained via a convolution of the in-medium 
nuclear matter g matrix with the target ground state mixed density. We have considered 
six different NN potential models in the construction of the corresponding omatrix, i.e. the 



II 

Paris potential jlQ] . the Nijmegen I, Nijmegen II and Reid 93 p otentials [l7j the Argonne 



AV18 potential [18| and the charge- dependent Bonn potential 19|. The g matrix is calculated 
in momentum space fully off shell, evaluated at eight different densities and nearly twenty 
momentum pairs, for all allowed NN states up to J = 7, in a square mesh of relative momenta 
up to 15 fm^^. The Fermi motion integrals involving the g matrix are made following 
Ref. Tj], where relativistic kinematics at the NN level are included in the folding integrals, 
supplemented by a proper treatment of the deuteron pole contribution 2^. For simplicity, 
and in order to keep our discussion focused within a single scheme, we have restricted the 
use of the NN charge-dependent potential models to their NN = components. 

The results reported in this study are based on semi-phenomenological densities. Proton 



densities were obtained from parametrized charge densities 2l|, unfolding the electromag- 



netic proton size to obtain the corresponding point density. Neutron densities were obtained 
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by adding to the proton density a harmonic oscillator orbital of the form ip ^ g-'^n^'^/s^ 
In the case of ^°Zr, a 3pG charge density was used, with the neutron excess characterized 
by a„=0.472 fm~^. This construction yields a proton (neutron) root-mean-square radius 
(Rrms) of 4.198 fm (4.363 fm). A similar construction was used for the ^'^Ca densities, with 
a 3pF charge density for the closed shells with the neutron excess characterized by a„=0.524 
fm^^. The resulting proton (neutron) Rrms is 3.374 fm (3.580 fm). In the case of ^^C, 
the proton density was represented by a modified harmonic oscillator. The neutron excess 
was characterized with a„=0.510 fm~^. The proton (neutron) Rrms is in this case 2.427 
fm (2.611 fm). The scattering results obtained with these densities were consistent with a 
variety of choices for the neutron-excess wavef unctions. 



A. Results for C{p,n) 



In Fig. [T] we present the measured 22j and calculated differential cross-section for the 
quasielastic ^^C{p,n) reaction following the (/-matrix approach described in the previous 



section. The excitation energy to the IAS, E^, is 2.31 MeV 23|. 

In this figure, and also in Figs. [1], [2] and IMTl we show results for the Paris (solid 
curves), CD-Bonn (dotted curves), Nijmegen I (short-dashed curves), Nijmegen II (long- 
dashed curves), Argonne AV18 (dot-dashed curves), and Reid-93 (dot-dot-dashed curves) 
NN potential models. Although in some cases the overlapping of the curves prevents their 
identification, this comparison allows one to visualize a moderate sensitivity of the calculated 
quantities to the choice of the A'^A'^ model. Such is the case of Fig. [H where the description 
of the differential cross section for scattering angles less than 25 degrees is reasonably well 
described by all NN potential models considered. A departure from the data is observed at 
larger angles, suggesting missing effects. 

Another quantity of particular interest is the so-called Fermi cross section {(Jf), that is 
to say, the forward-angle (zero-degree) differential cross section of the quasielastic {p,n) re- 
action. In Fig. [2] we present the measured and calculated Fermi cross sections for ^^C{p,n). 
The data are from Refs. [25[, and the various curves correspond to the six NN mod- 
els considered. A note of caution should be kept in mind before interpreting this figure. 
Although the (/-matrix full-folding optical model calculations account for relativistic kine- 
matics, all the NN potential models have been developed for energies below pion production 
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threshold. Thus, the results in Fig. [2] above 350 MeV represent extrapolations of the NN 
models. 

Below 150 MeV all NN models provide reasonably good descriptions of ap, albeit all of 
them underestimate the 80-MeV datum. As the energy increases, beginning at ~150 MeV 
results using the Paris potential depart from the others, being followed by those using the 
CD-Bonn potential near 300 MeV. Above pion production threshold the CD-Bonn potential 
exhibits a distinctive uniform growth, being followed less strongly by the Nijmegen I and 
Reid-93 models. Over the whole energy range both the Argonne AV18 and Nijmegen II 
potentials yield results in closest agreement with the data. In contrast, the Paris potential 
exhibits qualitative differences from the other models, particularly the depth of its minimum 
near 300 MeV. 

In a broad sense, the role of medium effects have been subject of significant interest in 
the description of NA processes. In Fig. [3] we contrast results for the Fermi cross section as 
function of the proton incident energy using the g matrix (solid curve) and t matrix (dashed 
curve). Here we also include results based on the off-shell tp approximation (dotted curve), 
a limiting case of the full-folding approach using the free t matrix and suppressing Fermi 
motion effects in the NN interaction. In all cases the Argonne AV18 potential is used. We 
observe that the differences between the g- and t-matrix approaches become more sizable 
at the lower energies. In general, the difference between the g— and t— matrix approaches 
above 150 MeV does not exceed 0.2 mb/sr, a useful upper bound of medium effects in 
{p,n) charge-exchange processes. At energies below 150 MeV the Fermi cross section exhibit 
sizable differences between both approaches. 



B. Results for ^^Ca(p,n) and ^'^^Zr{p,n) 



90 r 



In addition to ^^C{p,n), we have also performed calculations for quasielastic nucleon 



charge-exchange scattering on ^^Ca {E^ = 6.67 MeV |23j) and ^°Zr [E^ = 5.1 MeV |26j) 



igs. m and [5] we present results for '^^Ca{p,n) 



using the ^f-matrix full-folding approach. In 
at 135 and 160 MeV (data taken from Ref. 
we show results for ^^ZT{p,n) at 120 (data taken from Ref. 26|) and 160 (data taken from 



27|). respectively. Similarly, in Figs. [6] and [7] 



Ref. 



)MeV. The curve patterns follow the same convention as in Fig. [TJ 
In contrast to the ^^C(2),n) case, the Fermi cross sections for these four applications are 
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significantly underestimated relative to the data. The calculated ap for ^^Ca at 160 MeV 
averages 2.6(ib0.1) mb/sr, considering the six NN potential models. At the same energy, 
but in the case of ^°Zr, the average is 1.7(±0.1) mb/sr. The reported measurements for ^^Ca 
and ^°Zr are 5.0±0.3 mb/sr and 3.4±0.4 mb/sr, respectively, roughly twice as large as the 
calculated values. The measured-to-calculated ratio of ap at the lower energies increases to 
2.3 (^^Ca at 135 MeV) and 3.0 (^^Zr at 120 MeV). 

The results shown in Figs. HUT] correspond to fully consistent coupled-channel results 
involving explicitly the diagonal {Upp and f/„„) and off-diagonal {Ux) optical potentials based 
on the full-folding approach. Thus, they illustrate the degree of consistency among the 
three components of the potential within the theory. This is an important point to keep 
in mind, since it is possible to improve the description of the data with the use of the 
DWBA, supplemented by suitable phenomenological potentials for the nn channel. Indeed, 
in the upper frame of Fig. [8] we compare the differential cross section for proton scattering 
from ^°Zr at 121 MeV. The elastic scattering data are from Ref. 291], the solid curve 
represents the full-folding OMP results and the sparse crosses a standard Woods-Saxon 
(WS) parametrization. Although the full-folding results are in superior agreement with 
the data over the whole angular range, the quality of both descriptions at angles below 
25 degrees is comparable. A comparison at the hadronic level alone (i.e. suppressing the 
Coulomb term) is presented in the lower frame, where we show the differential cross section 
for uncharged nucleon scattering using the full-folding Upp (solid curve), f/„„ (dashed curve) 
and the local WS parametrization (sparse crosses). The similarity between the Upp and Unn 
is striking, suggesting their equivalence. The WS result is slightly larger than the other two 
cases. 

In a DWBA application it is customary to use, for the scattered neutron wave, Unn '■= Upp. 
In Fig. [9] we compare the consistent coupled-channel full-folding OMP results (solid curve) 
with the corresponding DWBA results (dashed curves). In both cases we use the same Ux 
based on the full-folding OMP. However, for the DWBA application we use the parametric 
WS potential already described. The significant increase of the differential cross section over 
whole the angular range is clear. 

In Fig. [To] we present the scattering amplitude (/) as function of the angular momentum 
/. The upper and lower frame are for the real and imaginary components, respectively. 
The filled (open) circles corresponds to exact coupled-channel results extracted from Eqs. 
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(137|) and (138!) . with the sohd curves representing the sum of both contributions. The stars 
(*) represent results based on DWBA using the full-folding x and nn couplings (DWBA- 
CC). The dashed curves correspond to DWBA calculations where the nn term has been 
replaced by a local WS parametrization (DWBA-WS). Clearly the imaginary component of 
the amplitude (lower frame) exhibits minor differences among the exact coupled-channel, 
DWBA-CC and DWBA-WS results. Additionally, we observe that the total amplitude is 
the result of a mutual cancellation of the Ux and Unn contributions. In the upper frame for 
the real components, however, we notice major differences between DWBA-WS and either of 
the two coupled-channel applications. Indeed, the real part of of the DWBA-WS amplitude 
is nearly twice as large as the other two. This difference is exclusively due to the nature 
of the diagonal coupling, extracted phenomenologically from other sources. In that sense, 
a phenomenological improvement of the data is achieved at the expense of a consistent 
two-channel description of the charge-exchange process. 

The difficulties discussed above have been explored in the context of the neutron-to-proton 
relative distribution. We have investigated the discrepancy between the consistent coupled- 
channel (yf-matrix OMP results and the measured quantities by studying the sensitivity of the 
scattering observables to variations in the neutron density |30|. Preliminary results along 
this line suggest the presence of a neutron halo in the nuclear periphery, that is to say, 
pronounced ratio of the neutron to proton density away from the nuclear surface. A more 
thorough investigation of this fact is a natural extension of the present work, and should 
not be limited to the radial matter distribution but include other features of the mixed 
density as well. An additional element in this discussion is the charge-dependence of the 
NN interaction, which in this work, has only been included in the = NN sector. 

C. Comparison with other results 

Various microscopic studies on {p,n) charge-exchange reactions have been presented in the 
past. The most recent one by Bauge et al. [31| consists of a Lane-consistent semi-microscopic 
optical model approach where an energy- dependent potential depth normalization factor is 
applied to the isovector components. These normalization factors are then used consistently 
in elastic and quasielastic scattering. Their quasielastic charge-exchange applications be- 
tween 100 and 200 MeV exhibit less pronounced discrepancies with the data relative to the 
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ones presented in this work. A microscopic extraction of these scahng factors remains to be 
seen. 

A relativistic description of {p,n) quasielastic reactions has also been provided by Clark 



et al. 



12] . In their case the relativistic impulse approximation is used, where Lorentz vector 



and scalar, isovector and isoscalar NA optical potentials are obtained from invariant NN 
amplitudes. Applications at 160 MeV for ^^Zi (p,n) lead to a Fermi cross section weaker 
by a factor of 2 relative to the data. Their interpretation is that the diagonal optical 
potential presents too much absorption. A sizable improvement in the description of the 
differential cross section is obtained with the use of phenomenological (diagonal) potentials 
in conjunction with DWBA. 

Although both approaches provide means by which to improve the description of the 
data, they also emphasize the limited understanding of this very simple process. 



V. SUMMARY AND CONCLUSIONS 



We have studied quasielastic {p,n) charge-exchange reactions to the IAS at proton en- 
ergies between 80 and 800 MeV. The study is based on an extension of the microscopic 
full-folding optical model approach for NA elastic scattering to processes where spin and 
isospin excitations become allowed. Thus, we have developed the necessary spin-isospin re- 
coupling in order to identify the off diagonal terms responsible for the exchange mechanism. 
In this particular work we focus on isospin excitations only. The resulting couplings, in the 
form of non-local optical potentials, are then applied to quasielastic {p,n) reactions on ^^C, 
^'^Ca and ^°Zr at intermediate proton energies. Emphasis is given to the full consistency 
of the microscopic approach, without resorting to the use of adjusting parameters of the 
model in order to describe the data. Indeed, it has been a primary focus of this work to 
disclose its limitations in its most complete current form, i.e. with an accurate account 
of the Fermi motion (off shell), and implied non localities in the optical potentials. Exact 
coupled-channel calculations, within numerical accuracy, were performed to obtain the scat- 
tering observables. Comparisons with the data show a persistent deficiency of the theory to 
account for them in a consistent fashion. This is particularly so in the case of the heavier 
targets, where the zero-degree Fermi cross sections are substantially underestimated. We 
observe, however, that the use of phenomenological representations of the optical potential 
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for the elastic channel yields sizable changes of the Fermi cross sections. However, this 
phenomenological correction does not allow a clear identification of the missing microscopic 
effects needed to account for the data. 

The above mentioned limitations become less dramatic in the case of {p,n) charge- 
exchange reactions on ^^C. In this case the full-folding optical model approach is able to 
explain reasonably well the differential cross section at 120 MeV. Additionally, depending 
on the NN potential model, the Fermi cross section is described with varying degrees of 
success in the energy range 80-800 MeV. These applications represent extrapolations of the 
of the underlying A'^A'^ potential model and favor the Nijmegen II and Argonne AV18 mod- 
els. A more definitive assessment of such results would require, however, the account for the 
above-threshold inelasticities at the NN level, an issue to be considered in the future. 

Though this work represents the most complete microscopic calculations to date, account- 
ing for full Fermi motion effects and implied non localities, with an accurate treatment of 
the resulting coupled-channel reaction equations, we conclude that the model is still unable 
to satisfactorily describe the data. The discrepancies are not minor, especially in the cases 
of {p,n) reactions on ^^Ca and ^°Zr studied here, where the zero-degree Fermi cross section is 
underestimated by factors of 2 or 3. Although the data can be described by using distorted 
waves based on phenomenological optical models, or by adjusting other parameters of the 
model, we have preferred to maintain the formal structure of the theory to stress the limita- 
tions of the microscopic model in its current form. There are other venues for the inclusion of 
additional effects not discussed here. Among the most immediate ones we mention the role 
of asymmetric nuclear matter in the NN effective interaction and the neutron distribution 
in the target ground state 30|] . Work toward inclusion of isospin asymmetry in the effective 
interaction is under way. 
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FIG. 2: The zero-degree Fermi cross section based on ^r-matrix full-folding optical model potentials 
using different NN potential models. For reference to the data sec text. The dashed line at 350 
MeV indicates the energy from which the calculations should be taken as extrapolations of all AW 
models. 
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FIG. 3: The zero-degree Fermi cross section based on full-folding OMP using the g matrix (solid 
curve) and t matrix (dashed curve). Off-shell tp results are represented with the dotted curve. All 
three results use the Argonne AV18 NN potential. 




23 



10 





:q 10 



^3 



10" 



_ 1 1 1 1 1 1 1 


1 ' : 


5 Ca(p,n) 


^'^J^ 160 MeV 


r \ ^ 








1.1,1. 


\ - 







10 15 
e [ deg ] 

cm. ° ^ 



20 25 



FIG. 5: Measured [27l] and calculated differential cross-section for ^^Ca(p, n) at 160 MeV. The 
curve patterns follow the same convention as in Fig. [2l 
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FIG. 6: Measured [2^ and calculated differential cross-section for ^^Zt{p, n) at 120 MeV. The curve 
patterns follow the same convention as in Fig. [5J 
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FIG. 8: Upper frame: Measured 
scattering from ^^Zr at 121 MeV. The solid curve corresponds to ^r-matrix full-folding OMP results, 
whereas the crosses correspond to a local WS parametrization. Lower frame: The calculated 
uncharged nucleon differential cross-section based on Upp (solid curve), (dotted curve) and the 
WS parametrization (sparse crosses). 
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FIG. 9: Measured 2a] and calculated differential cross section for ^^ZT{p,n) charge-exchange re- 
action based on full coupled-channel (solid curve) and DWBA using the parametric WS potential 
(dashed curve). 
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FIG. 10: The real (upper frame) and imaginary (lower frame) component of the scattering am- 
plitude / as function of the orbital angular momentum. The solid curves, * symbols and dashed 
curves represent the exact CC, DWBA-CC and DWBA-WS results, respectively. The filled and 
empty circles represent the Ux and [/„.„ contributions to the total amplitude [c.f. Eq. ([37|) ] . 
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